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ABSTRACT 


Rate Constants for the reaction of atomic chlorine 

with methane have been measured from 200-500K using the 

flash photolysis-resonance fluorescence technique,. When 

the results from fourteen equally spaced experimental 

determinations are plotted •' n Arrhenius form a definite 

curvature is noted « The results are best represented by 

16.34 _ ig 

the least squares equation, k = (5.44 + 4.08) x 10 

^(2,50 + 0.21) exp ( -608+62/T) cm^ molecule ^ s ^ „ The 

results are compared to previous work and are theoretically 


discussed 



INTRODUCTION 


Rate constants for the reaction of chlorine atoms with 

hydrogenous species are of great current interest due to the 

potential importance of such reactions in determining the 

. 1 

rate of ozone depletion in the earth’s stratosphere . Thus, 
termination of Cl atom chains in the stratosphere can occur 
via 

Cl + RH •+ HC l + R (1) 

which would result, at least temporarily, in the destruction 

of Cl atoms and interruption of the Cl-ClO chain which 

destroys ozone. The reaction of Cl atoms with CE 4 is of 

particular importance in this respect because of its 

abundance in the stratosphere^. 

Until very recently, only three direct studies had 

been made of the rate constant for the reaction 

Cl + CH 4 -» HC l + CH 3 ( 2 ) 

Davis, Braun and Bass 2 reported a value of k 2 - (1* 5+0.1) 

—13 3 "1 “1 

x 10 cm molecule s at 298K using the flash photolysis- 

3 

resonance fluorescence (FP-RF) technique. Clyne and Walker 

measured k g using a discharge flow - mass spectrometric 

(DF-MS) technique and obtained k^ = (5,08+0.53) x 

3 -1 -1 

exp (-17 90 +40 /T) cm molecule s from 300-686K, 


-11 

10 

while 


1 


4 

Poulet, LeBras and Combourieu , using a similar technique, 
reported kg=(l . 84+0 ,14) x 10 11 exp (-1400+100/T) cm 3 
molecule 1 s 1 from 295-490K. The room temperature values 
of k 2 from all three studies are in reasonable agreement, 
but the two temperature dependent studies clearly do not 
agree. 

Very recently several further studies have been performed. 

S -13 3 “X 

Leu and DeMore obtained k = (1.2+0. 3) x 10 cm molecule 

^ * 

-1 

s at 298K using the DF-MS technique, while the temperature 

dependence of kg has been further studied in several iabora- 

tories. Watson, Machado, Fischer and Davis, using 

_i o 

the FP-RF technique, obtained kg=(7, 94+0 .70) x 10 
exp ( -1260+3 5/T) cm. 3 molecule -1 s” 1 from 218-401K. The 
possibility of non-Arrhenius beha 'iour of k is clear from 
the above studies, in that A and fi are both higher at 

Ct £ 

higher temperatures, with reasonable agreement on the 
value of kg at room temperature. The possibility of a 

7 

curved Arrhenius plot was suggested by Manning and Kurylo 

—12 

who obtained kg=(7 . 93+1 . 53) x 10 exp (-1273+51/T) using 
FP-RF from 218-322K. However, they noted that their highest 
and lowest temperature points lay above the Arrhenius line. 


2 


while the mid -range points lay below the line, Their 

results from 277-322K could be fitted by the Arrhenius 

exjjression k 2 - (2,77 + 1,28) x 10 ^ exp (-1647 + 138/T) 

3 -1 -1 

cm molecule s which was close to that derived from 

8 

the results of Zahniser and Kaufman from 300-504K, The 
latter workers measured k^ from 200-504K using the dis- 
charge flow -resonance fluorescence (DF-RF) technique and 
observed non-linear Arrhenius behaviour. 

Much of the apparent discrepancies between the various 
Arrhenius equations reported for k g might be explained if 
the non-linear Arrhenius behavior obtained by Zahniser 
and Kaufman were confirmed using an independent technique. 
To this end, reaction (2) has been studied here at short 
temperature intervals in the range 200-500K using the 
FP-RF technique. 


3 


EXPERIMENTAL 


The FP-RF apparatus used in this study has been des- 

9 

cribed in detail previously. Consequently, only those 
aspects of the apparatus and procedures which are specific 
to the study of this reaction will be emphasized here. In 
all of the reported experiments, chlorine atoms were pro- 
duced by the flash photolysis of phosgene. 10 Except for 
one series of experiments at room temperature using a LiF 
filter, a sapphire filter was used to restrict photolyzing 
wavelengths to I ^ 143 nm so that photolysis of methane in 
the flash was negligible. Several preliminary experiments 
at room temperature using a suprasil filter with carbon 
tetrachloride as the atom source gave results identical to 
those reported here. 

Chlorine atom resonance radiation was produced by 
flowing a mixture of C-t^ i* 1 Se through the microwave 
discharge resonance lamp. The lamp output was observed 
with a McPherson Model 218 vacuum monochromator for a variety 
of compositions and pressure. Optimum conditions were 
established with 0.1% in He at a pressure of approximately 
0.5 torr. The multiple! structure showed only moderate 
reversal under these operating conditions. 11 Resonantly 
scattered photons were viewed at right angles through a 


4 
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Bar filter without wavelength resolution, and the signal 
2 

was assumed to be ' inearly proportional to the atom 
12 

concentration,, The 3-component reaction mixtures of 
methane, phosgene source and argon diluent flowed through 
the cell at a rate sufficient to replenish the reaction 
mixture between successive flashes. Preliminary experiments 
showed, that pseudo-first -order rate constants, obtained 
from exponential decay of the resonance fluorescence, were 
slightly higher under either static or very slow flow 
conditions, indicating that secondary complications could 
contribute at insufficiently fast flow rates; 

The pseudo-first -order rate constants obtained in the 

presence of a large excess of methane were composite since, 

in addition to reaction with CH^ and any impurities, Ct atoms 

were lost by diffusion out of the reaction viewing zone. 

The diffusion correction, k^, was determined independently 

9 

in the normal way by flash photolyzing mixtures of phosgene 
and argon at the various experimental temperatures and 
pressures employed to obtain the kinetic results. At very 
high incident flash intensities, i,e, at high flash energies 
and with the optical train from the flash lamp to the reaction 
cell freshly cleaned, k d was observed to depend upon the 
intensity of the flash. Such variations were easily allowed 
for in practice. The k^ correction- term under such conditions 



never exceeded 10% of the observed decay rate constant and 


was normally « 10%, For most series of experiments was 

independent of intensity and was again typically «10%, 
although it was -15% at the highest temperatures studied 
here , 

Argon (Matheson, 99,9995%) and helium (Airco, 99,9999%) 

were used without further purification. Chlorine (Matheson, 

99,5%) was further purified by fractional distillation at 

195K and phosgene (Matheson, 99%) was similarly purified 

by fractionation at 163K, Methane (Matheson, 99,97%) was 

used without further purification. Gas chromatographic 

analysis of the CH. showed it to contain 32ppm of C H 

4 ^ D 

as impurity, but no other hydrocarbon impurities were 
detectable. 


* 


6 



RESULTS 


Under the pseudo-first-order conditions employed here 
with [CH^] » [C-t], the decay of Ct atoms may be represented 


by 


tn [Ct] = -k . . t + -tn [Ct] (3) 

observed o 


where the observed pseudo-first -order rate constant is 
given by 

k - , “ k [CH. ] + k . (4) 

observed bi 4 d 

Plots of £n [C-t] vs t were linear and k , ,, and hence 

observed' 

k fai , was obtained from such plots using a linear least 

squares analysis. The results are presented in Table 1, 

and as shown there, k^ was independent of substantial 

variations in both total pressure and [CH^] . The essential 

correctness of equation (4) is indicated by the linearity 

of plots of k . , vs [CH.], as shown in Fig. 1. 

^ observed 4 ’ & 

Despite the very high [stable reactant] / [atom] ratios 

13,14 

accessible in this system, it has been shown previously 
that contributions from secondary reactions resulting in 
additional atom loss can become important at very high 
incident flash intensities i„e„, at relatively high atom 
and/or free radical concentrations. Consequently, the 
reaction was studied over a wide range of flash intensity, 
and hence a wide range of [C -ll, under all conditions. For 
most experimental conditions, the measured k^ was 


7 


independent of the accessible flash intensity range. However, 
for some series of experiments, in particular for those 
performed with a freshly cleaned optical train from flash 
lamp to reaction cell, the resulting was observed to 
increase slowly from its low flash energy limit to ~ , 20% 
above that limit at very high flash energies. This increase 
in was indicative of secondary complications and, in 
the few series of experiments where this occurred (at and 
around room temperature) , only those experiments which were 
clearly independent of a range of intensity are included in 
Table 1. 


In order to obtain k 2 from the values listed in 
Table 1, corrections to k^ must be applied to account 

7 

for the fast reaction of Ct with the 32ppm CgH^ impurity 

always present in the CH^ used here. Table 1 includes a 

listing of lc bi ( correoted )» i ;he latter being equal to kg. 

The correction decreases from ^10% at 200K to ~o„2% at 500K. 

Fig. 2 shows an Arrhenius plot of the mean kg obtained 

at each of the 14 temperatures. A linear least squares 

—12 

treatment yields the equation kg^Cll .0+1 .2) x 10 

3-1-1 

exp(-1350+30/T) cm molecule s . However, a definite 
curvature is apparent in Fig. 2. Thus the results may be 


8 


better represented by an equation of the form k~AT n exp(-B/T) 

and the least squares fit, shown as the dashed line in Fig. 2 

...... . . . ft . n - n — !9 (2 . 50+0 . 21) 

is given by the equation kg- (5.44 -4.08) x 10 T — 

exp (-60 8+6 2 /T) cm 3 molecule 1 s \ 


9 


I 


1 
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DISCUSSION 

The apparent uncertainty in the Arrhenius parameters 

for k 0 is evident in the very recent review compiled by 

15 -11 

Watson. The reported A 2 values range from 5.08 x 10 

-12 3 —1 —1 

to 7.0 x 10 cm molecule s while the Eg values range 

-1 

from 3.56 to 2.50 kcal mol . It is possible that non- 
linear Arrhenius behaviour may* at least in part, be 
responsible for these discrepancies. This has been suggested 

7 

by Manning and Kurylo and is reinforced by the results ' 

obtained here. Thus ft the present results are analyzed 

temperature ranges equivalent to those used in other 

studies, as shown in Table 2, the agreement with other 

studies using resonance fluorescence detection is excellent. 

4 

The agreement with the DF-MS study of Poulet et al also 

is good. However this may be fortuitous, given the poor 

precision of the data from reference (4) . The other DF-MS 

3 

study of Clyne and Walker yields Arrhenius parameters 
significantly higher than the resonance fluorescence 
results. Nevertheless, it is clear that much of the apparent 
conflict in Arrhenius parameters of different studies may 
be explained by non-linear Arrhenius behavior, with the 
resulting A„ and E depending on the T range. 


10 


uhBIUTY OF THE 


The discrepancy between the DF-MS results, obtained 

under Ct-rich conditions, and those obtained using resonance 

fluorescence detection under CH^-rich conditions, is presented 

graphically in Fig. 3. The DF-MS results, represented by 

filled symbols, seem generally ~ 20% higher than the open 

symbol? of the resonance fluorescence data and in general 

the two sets of results do not. quite overlap within the 

claimed uncertainties. Nevertheless it is clear from Fig. 3 

that the rate constant is known with a great deal more 

precision than most gas kinetic rate constants. This is 

particularly so between 200-300K, since all of the resonance 

fluorescence results, both DF-RF and FP-RF, are in remarkably 

good agreement. The present results are essentially identical 

7 

to those of Manning and Kurylo and lie between the slightly 

0 

higher values of Watson et al and the slightly lower values 

8 

of Zahniser and Kaufman. 

However, as indicated in Fig. 2 and Table 2, the present 

results over the complete 200-500K range seem best represented 

+16.34 —19 

by a 3-parameter fit given by k =(5.44- 4.08) x 10 
T (2.50+0.21) exp (-608+62/T) cm molecule s . All of the 

resonance fluorescence results from references (2) , (6) , 

(7), and (8), together with the present results, can be 

+88.32 _17 

fitted to the similar expression lc 2 =(2.08«2.03 ) x 10 


11 



T (1.96^0. 56) exp (-766+168/T) cm^ molecule ^ s ^ for 200- 
504K, The difference between these two equati^”? is 
largely a result of the uneven distribution of the combined 
data over this temperature range, in particular the paucity 
of data at T>407K. 

A theoretical description of reaction (2) has already 

been reported with the BEBO method in order to explain 

measured kinetic isotope effects, ^ and these calculations 

included one dimensional tunneling corrections. We have 

extended the BEBO calculation by including triplet repulsion 

17 

and using separate values for bond indices. Also one 

dimensional Eckart tunneling factors have been included 

18 

after the method of Truhlar and Kuppermann in contrast 

17 19 

to that of Johnston, ’ The two methyl wags in the 

activated complex are assumed to be degenerate, and its 

frequency is taken as a parameter in fitting the data. For 

a value of 1140 cm \ the computed values are well 

“19 2 5 

represented by k = 8.4 x 10 T /tU exp(-748/T) over the 

temperature range 200-500K, 

The agreement between calculated and experimental 

results is excellent but may be fortuitous given the 

approximate nature of BEBO and the uncertainties surrounding 

19 

both the application of one dimensional tunneling theory 


12 


r 


T 


T 


and the ambiguities of activated complex vibration frequency 
assignments * These inadequacies are further reflected in 
isotope effect calculations,, Consistent extension of the 
model to Cl + CD^ predicts (k H /kp)300K ~ M , whereas the 
measured value is 11.0* lo Thus, the experimentally 
observed non-Arrhenius behavior cannot be used to imply 
the correctness of the model and, in particular, the 
presence of tunneling* On the other hand the model as 
described is consistent with experiment* We point out for 
the present case that the greatest ambiguities arise from 
tunneling corrections and vibration frequency assignments 
and not from the method of obtaining the potential energy 
of interaction (BEBO) , Even with more sophisticated methods 
(eg* LEPS) , the above mentioned ambiguities will still be 
present and will make the theoi’etical results somewhat 
arbitrary* 
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Table 1. Rate data for the reaction Ct + CH4 -* KCt + 
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Table 2. Linear least squares Arrhenius parameters for k g 



T range 

A 2 a 

E 2 a 

Reference 

K 

-12 3 _ _ -1 -1 

10 cm molecule s 

kcal mole -1 



218-322 

7 . 93+1 o 53 

2.53+0.10 

7 

22C-318 

7 . 83+0 . 67 

2 o 54+0 o 04 

this work 

218-401 

7.94+0.70 

2 ,> 50+0 1 07 

6 

220-404 

9.91+0.67 

2.66+0.04 

this work 

200-504 

13.6+2,0 

2.87+0.08 

8 

200-500 

11 o 0+1 . 2 

2.69+0.06 

this work 

296-504 

22 .5+3.1 

3.23+0.08 

8 

299-500 

18.4+2.8 

3.07+0.10 

this work 

296-490 

19.2+5.8 b 

2.86+0.19 

4 

300-484 

42.4+14.1° 

3.42+0.22° 

3 


a error limit is the standard deviation . 

b recalculated from the original data; the author’s quoted 
A-factor and error limit for A appear to be wrongly 
calculated in ref. 4„ 

c although this study was from 300-686K, these parameters 
apply only to the data from 300-484K. 




FIGURE CAPTIONS 


Fig. 1. Pseudo -first -order rate constant k . . vs CCH„] for 

observed 4 

the conditions O: 276 K, “ 5 0 Torr; •: 299 K, 

]?total - 50 Torr, and displaced upwards by 100s '• 318 K, 

P total 5=1 50 Torr > and displaced upwards by 200s’" 1 . 

Fig. 2. Plot of log 10 k 2 (-log 10 (corrected) vsi The 

dashed line represents log.^ k 2 e valuated from the equation 
kg = 5. '44 x 10 1J * t 2 ' 50 exp(-608/T) cm^ molecule 1 s’" 1 . 

Fig. 3. Comparison of log 1Q k 2 va i ues > filled symbols represent 
DF-MS data while open symbols are DF-RF and FP-RF data, 
as follows: •, ref. 3; ▼, ref. 5; A, ref. 4 ; o, ref. 
v ; ref. 2; A, ref, 6; 0* , ref. 7 j 0, this work. 
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